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In MHD symmetric systems the equilibrium physical quantities are dependent on two variables 
only. In this cases it is possible to find a magnetic surface function that has the same symmetry. 
Under the assumption that the metric determinant is also independent of a third, ignorable coordi- 
nate, a general MHD equilibrium equation in curvilinear coordinates is deduced. This equation is 
specially useful when non-orthogonal generalized coordinates are used. 
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I. INTRODUCTION 



II. CURVILINEAR COORDINATES 



Ideal magnetohydrodynamics (MHD) is the most basic 
single-fluid model for determining the macroscopic equi- 
librium and stability properties of a plasma [TJ [5] . There 
is enough evidence that this model describes how mag- 
netic, inertial and pressure forces interact within a per- 
fectly conducting plasma [3J H] . 

The ideal MHD is used in magnetic fusion to describe 
static equilibria and to infer convenient magnetic geome- 
tries for confinement. In symmetric plasma systems the 
field lines lie on a set of closed nested toroidal magnetic 
surfaces. These surfaces can be determined by solving 
the MHD equilibrium equation E] . 

Equilibrium equations for different symmetric plasma 
systems have appeared in the literature [2]- In partic- 
ular, the well known Grad-Shafranov [7] equation is 
written in terms of orthogonal coordinates and describes 
axisymmetric toroidal equilibrium. 

In this article a general MHD symmetric system is de- 
scribed by a magnetic surface with the same symmetry 
that the considered equilibrium. This function is ob- 
tained from a general equilibrium equation in curvilinear 
coordinates. With this equation, once the curvilinear co- 
ordinates has been chosen, the equilibrium equation in 
any geometry can be derived. This is especially useful 
when non-orthogonal generalized coordinates are used. 

The general curvilinear coordinates used in this ar- 
ticle are defined in Section I. A symmetric transversal 
magnetic flux VP and a function / that determines the 
transversal electric current are introduced in Sections II 
and III. In Section IV, a general equilibrium equation re- 
lating the equilibrium pressure to the surface functions 
^ and / is derived. Finally systems well-known in the 
literature are considered as examples and the equilib- 
rium equations in terms of conventional toroidal coor- 
dinates (Appendix A), helical coordinates (Appendix B) 
and natural coordinates (Appendix C) are derived from 
the general equation obtained in this article. 
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In symmetric plasma confinement systems, all the equi- 
librium functions having a physical meaning are depen- 
dent on two variables only. Curvilinear coordinates are 
named u\,U2 and u%. The surfaces Ui — Ci, where Ci is a 
constant, are coordinate surfaces. A coordinate curve itfc 
is a curve along which Ui, uj (i ^ j ^ k) are constants. 

The coordinates u\ and U2 are chosen in order to have 
the magnetic axis of the system coincident with a coordi- 
nate curve 1*3 and U2 is a transversal coordinate, 1*3 will 
be an ignorable coordinate; longitudinal directions are 
given by coordinate curves u\. In plasma confinement 
problems we have, usually, periodicity in it 2 and W3. The 
following periodicity 



u 3 = L(ui,u 2 ) 



(1) 



is assumed. 

An attempt is made use of the notations most familiar 
in the literature. 

The covariant basis vectors are given by 



Or 

dui 



(2) 



where is tangent to the Ui curve and the contravariant 
basis vectors are defined by 



e l = Vit, 



(3) 



where e l is normal to the Ui surface. U\,U2 and U3 are 
taken in order to satisfy 



e 4 = y/ge 1 x e k 



(4) 



for any cyclic permutation (i,j,k), g is the determinant 
of the covariant metric tensor g^ j. 



III. TRANSVERSAL MAGNETIC FLUX 

Define M{ui, u 2 ) as the magnetic flux through a coor- 
dinate surface 112 which extends from the magnetic axis 
to a coordinate curve U3 and limited by < 113 < L. On 
the magnetic axis U\ = a and 



B = B 3 e 3 



(5) 
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Then: 



dv! Y 



from which it follows 
<9* 



1 



dui L J 
Taking account of the equation 
V • B = 



fgB~du 3 



\/gB 2 du 3 



(6) 



(7) 



(8) 



and assuming B 1 = on the axis we can derive from 
expression (JTj]): 



du 2 



^fg~B 1 du 3 



(9) 



If ^fgB 1 and *JgB 2 are independent of u 3 , we find an 
expression for B in terms of \l/: 



e 3 e 3 
S = — x W + B 3 — 



(10) 



533 533 

^ = constant represents a magnetic surface because 

V*-B = (11) 



as can be seen using (10). 

The magnetic flux can also be expressed in terms of 
the vector potential A using the Stokes theorem: 



A 3 du 3 , 



(12) 



A 3 is assumed to be zero on the axis. In symmetric sys- 
tems: 



$(ui,u 2 ) = -A 3 (ui,u 2 ). 



(13) 



The only restriction to the gauge of A is in order to keep 
the same symmetry as the physical quantities. 



IV. CURRENT DENSITY 

The current density satisfies the equations 

VxB = fi J and V • J = 0. (14) 

On the magnetic axis J 1 = 0. These equations are similar 
to the equations for B: 



Vx4 = BandV-B = 



(15) 



and also B 1 — on the axis. 

Similar considerations must yield similar results. Let 
us define a function: 



Vol = J B 3 du 3 , 



(16) 



which is an expression similar to (12). In a symmetric 
case it would be: 



Vol = -B 3 (ui,u 2 ). 

In this case on the axis B 3 is not zero. 
The transversal current is given by: 



L(I - I a 



yfgj 2 du\du 3 . 



(17) 



(18) 



This expression can be compared to (f6j. Thus consid- 
ering the symmetry argument, the current density is ex- 
pressed in terms of I as: 



T e 3 e 3 
J = — x V/ + J 3 — . 

533 533 



(19) 



an expression similar to ( 10 1 



V. PRESSURE EQUILIBRIUM EQUATION 

The MHD equilibrium theory, with scalar pressure P, 
considers the equations: 



VP = J x B 



and 



Vx6 = [m J, 
L P and 9 satisfies the relations: 



B ■ V* 
J ■ VP 



B ■ VP 

J V/ = 



= 
0. 



(20) 

(21) 

(22) 
(23) 



what means that /, P and '5 are surface quantities. 
Using (20) together with (17) and (19) a relation be- 



tween these surface quantities is found: 



VP 



933 



B3 

933 



-VI, 



(24) 



J 3 can be taken from (|21| and P 3 
pi| and u! 

(A*^)V«- 



-Hoi. Substituting 



them in ([24j) and using ( 10 ) we find the final expression: 
-M0.933VP - fJ, 2 IVI 



where 



533 


' d 


( 523 \ 


d 


( 513 


Vs 




V533/ 


du 2 


V533 



V*, (25) 



533 



y/9 g 33 
d V9 ( 12 3* 



d ^9 ( u a* 

9 8 Ul 



,89 
du 2 



du 2 g 33 



9 



dui 



9 



du 2 J 



(26) 



/ and P are functions of only. Therefore, whenever 
V*? 7^ 0, the expression can be simplified to a scalar 
equation 



A** = -n g 33 P' -ij, 2 II' 



933 


' d 


( 923 \ 


d 


( 513 


V9 




\933J 


du 2 


V533 



(27) 
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which corresponds to the Grad-Shafranov equation [2J. 
Here the prime indicates differentiation with respect to 
"J. The quantity ^g^P + Mo^ 2 /2 must be continuous 
through a surface where = 0. 

The components of the magnetic field are expressed in 
terms of ^> as: 



V9B 2 = 
and 



du 1 



(28) 



V9B 1 = - 



du 2 



The equation (26) can also be written as: 



A** = V 2 *-W-— ^ 

533 



or 



A** = 533 V 



■(-) 

V533 ) 



(29) 
(30) 



The general equilibrium equation can be specially useful 
when non-orthogonal generalized coordinates are used. 

In Appendix A toroidal pinches with axial symmetry 
are considered, and the equilibrium equation in conven- 
tional toroidal coordinates [H] is derived from equation 
(14). If the system presents a straight helical symmetry, 
the coordinates u\ = r, u 2 — — az and u 3 — z can be 
introduced, where a is the pitch of the helix. Using ( |25[ ) 
the equilibrium equation as found on the literature [2] 
is deduced straightforwardly. In these systems it is very 
likely to appear discontinuity surfaces where = 
(See Appendix B). The well-known equilibrium equation 
using flux coordinates [2] is obtained also very easily in 
Appendix C. 



VI. CONCLUSIONS 

In this article a general MHD symmetric system has 
been considered, and a magnetic surface function, with 
the same symmetry, has been introduced to describe it. 
An equilibrium equation satisfied by this function is de- 
duced in curvilinear coordinates. This equation is a gen- 
eralization of several particular MHD equilibrium equa- 
tions valid for the equilibria considered in the literature. 
Therefore, given the curvilinear coordinates, the equilib- 
rium equation for a magnetic surface function, in any 
geometry, can be derived. This is especially useful if 
nonorthogonal coordinates are used. As examples, in the 
appendices, equilibrium equations well known in the liter- 
ature are obtained from the mentioned general equation 
presented in this article. The procedure followed in this 
paper resembles the one used to deduce an equilibrium 
equation for incompresible irrotational steady fluid flow 



Appendix A: Toroidal pinch with axial symmetry 

The equilibrium equation using a conventional toroidal 
coordinate system [8] (see Figjl]). 

"1 = £ , u 2 = w , u 3 = ip 

is obtained in this appendix. The coordinates are defined 
by: 



i?Q sinh £ 
cosh£ — cosw 



R' Q sin uj 
cosh £ — cos uj 



were r, z and ip are the polar cylindrical coordinates and 
Ro is the major radius. If £ = £0 defines the toroidal 
surface, then cosh£ = Ro/b, R' = i?o\A — b 2 / R$. 




FIG. 1: Toroidal coordinate system. 
The contravariant basis is 



e 3 = Vip 



with 



and, 



he = hu = 



h 



R'n 



cosh £ — cos ui 



hip = /i£ sinh £ 
, e w and e v are unit vectors. The metric is given by: 
y/g = (e 1 ■ e 2 x e 3 ) -1 = h^h^ 

11 1 22 j.2 

9 = Jp,=9 ; 333 = V 



The equilibrium equation (27) becomes: 



sinh£ / 9 1 W d 1 9* 



9, 



ti£ \ <9£ h v d£ doj h v doj 



(Al) 
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If a function F is introduced as: 

* = [2(cosh£ - cosuj)]- 1/2 F (A2) 
we obtain the well known Grad-Shafranov equation [2J: 



d 2 F d 2 F , BF 1 

f 4 J R' 2 sinh 2 £ , , , 

X { [2(coshVcost)] 2 ^ P+ ^ / 



2R' 2 



2(cosh£ — cosw)] 3 / 2 
(A3) 



The magnetic held is derived using (27 1: 

<9* 



and 



h^h^pBu = 



<9w 



(A4) 



Appendix B: Helical system with a straight 
magnetic axis 



The coordinates are: 
u\ = r , u 2 - 



az — u , u 3 



where r, and 2 are the polar coordinates. Their con- 
travariant basis is: 



1 2 
e = e r ,e 



r 



where e ri eg and e z are the unit vectors in cylindrical 
coordinates, the covariant basis is: 

ei = e r , e 2 = re e , e 3 = e z + aree 

The metric is given by: 



9 



13 — 



1 

^+a 2 -a 
-a 1 



1 

r 2 ar 2 
ar 2 1 + a 2 r 2 



In this case the equilibrium equation ( 27 ) becomes [2] 



1 + err 



2„2 



9 



Id 2 * 



dr \ 1 + a 2 r 2 9r / r 5m 2 

1 + a 2 r 2 d dr 2 



(l + a 2 r 2 V P'- M 2 //' + Mo 7 



dr 1 + a 2 r 2 
(Bl) 



and the magnetic held components are derived using the 
following equations: 



<9* 

dr 



Ba - arB„ 



du 



-rB r 



(B2) 



Appendix C: Natural coordinates 

In natural systems [10] the coordinate u\ is a magnetic 
surface label, it is analogous to the minor radius of the 
torus, u 2 and 113 are poloidal and toroidal cyclic coordi- 
nates, ranging from to 27T. 

The physical variables are: 

Magnetic pressure P(u\) 

Poloidal flux 2wx(ui) = Jq 1 du[ f Q n y/gB 2 du s 
Poloidal current 27r(J — I axis) — Jq 1 du'i 

x Jo * V9 j2du 3 
Toroidal flux 2tt$ = f^ 1 du\ y/gB 3 du 2 

Toroidal current 2uJ = J Q 1 du' x J Q * y/gJ 3 du 2 

Volume inside a (2n) 2 V(ui) — f™ 1 du^ J 2ir du 2 

toroidal magnetic x f * ^fgduz 
surface 



The equilibrium equation ( 24 ) becomes 



P = x 1 



333 



333 



(CI) 



where the prime indicates derivation with respect to u%. 
Multiplying equation (CI) by y/gdu 2 du3 and integrating 
in a magnetic surface we get 

P'V = -x'J' - &I' (C2) 

The magnetic field and the current density are given by: 



333 



333 



and 



333 



333 



(C3) 



(C4) 



In terms of the contravariant components ( C4 1 becomes: 

(C5) 



e 2 + J 3 e 3 



Introducing the average value of J 3 yfg in a magnetic sur- 
face: 



1 

2^ 



J 3 y/gdu 2 



and writing 



expression (C5) becomes 



I' J' J 3 Jg 
J = —e 2 + — e 3 H ^-e 3 

•Jg \[g V9 
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Defining a function v by: 
dv 



where $ is the scalar potential in the absence of the 



plasma. This equation equated to ( C3 ) represents the 



du 2 



= J 3 Vg 



the current density can be written: 

J = V x (-Ie 3 + Je 2 - ve 1 ) 
and B comes naturally as: 



equations commonly used to determine the flux coordi- 
nates [HJ1ITT]. 

The only assumption taken is of the symmetry of the 
system. 



B 

Mo 



= -Ie 3 + Je 2 - ve 1 + V$ 
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